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Abstract. In this paper, we define a stochastic Sierpinski gasket, on the basis of which we construct a 
network called random Sierpinski network (RSN). We investigate analytically or numerically the statistical 
characteristics of RSN. The obtained results reveal that the properties of RSN is particularly rich, it is 
simultaneously scale-free, small-world, uncorrelated, modular, and maximal planar. All obtained analytical 
predictions are successfully contrasted with extensive numerical simulations. Our network representation 
method could be applied to study the complexity of some real systems in biological and information fields. 

PACS. 89. 75. He Networks and genealogical trees - 89.75.Fb Structures and organization in complex sys- 
tems - 05.10.-a Computational methods in statistical physics and nonlinear dynamics - 87. 23. Kg Dynamics 
of evolution 



1 Introduction 

' In the last few years, much attention has been paid to the 
study of complex networks as an interdisciplinary sub- 
ject pp. It is now established that network science is a 
powerful tool in the analysis of real-life complex systems 

■ by providing intuitive and useful representations for net- 
' worked systems. Many real-world natural and man-made 

systems have been examined from the perspective of com- 
plex network theory. Commonly cited examples include 
| the Internet [2], the World Wide Web [3], metabolic net- 

■ works [4], protein networks in the cell [5], co-author net- 
' works [6], sexual networks 7J, to name but a few. The 

empirical studies have uncovered the presence of several 
generic properties shared by a lot of real systems: power- 
law degree distribution [H], small- world effect including 
small average path length (APL) and high clustering coef- 
ficient [9] , and community (modular) structure |10j . These 
new discoveries have inspired researchers to develop a va- 
riety of techniques and models in an effort to understand 
or predict the behavior of real systems [Tj . It is still of cur- 
rent interest to reveal other different processes in real-life 
systems that may lead to above general characteristics. 

In our earlier paper, we have proposed a family of 
deterministic networks based on the well-known Sierpin- 
ski fractals [11]. These networks posses good topologi- 
cal properties observed in some real systems. However, 
their deterministic construction are not in line with the 
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randomness of many real-world systems. In this paper, 
we present a stochastic Sierpinski gasket, in relation to 
which a novel network, named random Sierpinski network 
(RSN), is constructed. The obtained network is a maximal 
planar graph, it displays the general topological features 
of real systems: heavy-tailed degree distribution, small- 
world effect, and modular structure. We also obtain the 
degree correlations of RSN. All theoretical predictions are 
successfully confirmed by numerical simulations. 



2 Brief introduction to Sierpinski network 

In our previous work [11] , motivated by the classic de- 
terministic fractal, Sierpinski gasket (or Sierpinski trian- 
gle) [T2"I[T3"] , we introduce a new type of graph, called Sier- 
pinski network. The well-known Sierpinski gasket, shown 
in figure [Ha), is constructed as follows [T3]: We start with 
an equilateral triangle, and denote this initial configura- 
tion by generation t = 0. Then in the first generation 
t = 1, the three sides of the equilateral triangle are bi- 
sected and the central triangle removed. This forms three 
copies of the original triangle, and the procedure is re- 
peated indefinitely for all the new copies. In the limit of 
infinite t generations, we obtain the famous Sierpinski gas- 
ket, whose Hausdorff dimension is df = In 3/ In 2 [14]. 

From this famous fractal we have defined the Sierpinski 
network [11] as illustrated in figure QJb). The translation 
from the fractal to network generation is quite straight- 
forward. Let each of the three sides of a removed triangle 
correspond to a node (vertex) of the network and make 
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Fig. 1. The first two stages of construction of the Sierpinski 
gasket (a) and its corresponding network (b). 

two nodes connected if the corresponding sides contact one 
another. For uniformity, the three sides of the initial equi- 
lateral triangle at step also correspond to three different 
nodes. The resultant Sierpinski network has a power-law 
degree distribution P{k) ~ fc~ 7 with 7=1 + j 2 ^, displays 
small-world effect, and is disassortative. 

3 Random Sierpinski network and its iterative 
algorithm 

In this section, we first construct a random Sierpinski gas- 
ket from the deterministic Sierpinski gasket (or Sierpinski 
triangle) [13] . Then we will establish a random Sierpinski 
network based on the proposed stochastic fractal. Analo- 
gous to the Sierpinski triangle, the random Sierpinski gas- 
ket also starts with an equilateral triangle. At step 1, we 
perform a bisection of the sides and remove the downward 
pointing triangle forming three small copies of the original 
triangle. Then in each of the subsequent generations, an 
equilateral triangle is chosen randomly, for which bisection 
and removal are performed to form three small copies of 
it. The sketch map for the random fractal is shown in the 
upper panel of figure [U From this fractal we can easily 
construct the random Sierpinski network with sides of the 
removed triangles mapped to nodes and contact to links 
between nodes. As in the construction of the deterministic 
version [11] . the three sides of the initial equilateral tri- 
angle at step are also mapped to three different nodes. 
Figure [2] (lower panel) shows a network derived from the 
random Sierpinski gasket. 

According to the construction of random Sierpinski 
network, we introduce a general iterative algorithm gener- 
ating the network. We denote the random Sierpinski net- 
work after t iterations by W(t), t > 0. Initially (t — 0), 
W(0) has three nodes forming a triangle. At step t = 1, 
we add three nodes into the original triangle. These three 
new nodes are connected to one another shaping a new 
triangle, and both ends of each edge of the new triangle 
are linked to a node of the original triangle. Thus we ob- 
tain W{1). For t > 1, W(t) is obtained from W(t - 1). 
For the convenience of description, we give the following 
definition: For each of the existing triangles in W[t — 1), if 
there is no nodes in its interior and among its three nodes 
there is only one youngest node (i.e., the other two are 
strictly elder than it), we call it an active triangle. At step 
t — 1, we select at random an existing active triangle and 




Fig. 2. (Color online) The sketch maps for the construction of 
random Sierpinski gasket (upper panel) and its corresponding 
network (lower panel). 




Fig. 3. (Color online) Iterative construction method for the 
network. 



replace it by the connected cluster on the right of figure [3] 
then W(t) is produced. 

Since at each time step the numbers of the nodes and 
edges increase by 3 and 9, respectively, we can easily know 
that at step t, the network consists of N t — 3t + 3 nodes 
and E t = 9t + 3 edges. Thus, the relation E t = 3N t — 6 
holds for all steps. In addition, according to the connec- 
tion rule, arbitrary two edges in the network never cross 
each other. Therefore, the considered network is a maxi- 
mal planar graph [15] , which is similar to its deterministic 
version [11] and some previously studied networks [16] , 
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4 Structural characteristics 

We now study the statistical properties of RSN, in terms 
of degree distribution, clustering coefficient, average path 
length, degree-degree correlations, and modularity. The 
analytical approaches are completely deferent from those 
applied to the deterministic Sierpinski network [llj . 



4.1 degree distribution 

Initially (t = 0), there is only one active triangle in the 
network. In the subsequent iterations, at each time step, 
three active triangles are created and one active triangle is 
deactivated simultaneously, so the total number of active 
triangles increases by 2. Then at time t, there are 2t + 1 
active triangles in RSN. Note that, for an arbitrary given 
node, when it is born, it has a degree of 4 and one active 
triangle containing itself; and in the following steps, each 
of its two new neighbors separately generates a new active 
triangle involving it, and one of its existing active triangles 
is deactivated at the same time. So, for a node with degree 
k, the number of active triangles containing it is Let 
Nk (t) denote the average number of nodes with degree k at 
time t. By the very construction of RSN, the rate equation 
that accounts for the evolution of N k (t) with time t is [T7] 



dN k (t) 
dt 



fc-4 



N k - 2 (t) 



k-2 



N k (t) 



2t+l 



36 



k,4- 



(1) 



The first term on the right-hand side (rhs) of Eq. (fT]) ac- 
counts for the process in which a node with k — 2 links is 
connected to two new nodes, leading to a gain in the num- 
ber of nodes with k links. Since there are N k -2{t) nodes of 
degree k — 2, such processes occur at a rate proportional to 
^=iVfe_2(£), while the factor 2t + 1 converts this rate into 
a normalized probability. A corresponding role is played 
by the second (loss) term on the rhs of Eq. ([1]). The last 
term on the rhs of Eq. (fT]) accounts for the continuous 
introduction of three new nodes with degree four. 

In the asymptotic limit Nk(t) = (3t + 3)P(k), where 
P(k) is the degree distribution. Substitute this relation 
into Eq. (fTJ) to lead to the following recursive equation for 
infinite t 



k > 4 
for k — 4 , 



giving 



P(k) = 



32 



k(k + 2)(k-2) 



(2) 



(3) 



In the limit of large k, P(k) ~ fc~ 3 , which has the same 
degree exponent as the BA model [5] and some hierarchi- 
cal lattice models [18]. In order to confirm the analytical 
prediction, we performed numerical simulations of the net- 
work plotted in figure [¥J which shows that the simulation 
result is well in agreement with the analytic one. 
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Fig. 4. (Color online) Log-log graph of the degree distribution 
for a network with order N = 300000. The squares denote 
the numerical results and the solid line shows the theoretical 
predication given by Eq. Q. 



4.2 clustering coefficient 



By definition, the clustering coefficient [9j Ct of node i is 
defined as the ratio between the number of edges e; that 
actually exist among the ki neighbors of node i and its 
maximum possible value, ki(ki — l)/2, i.e., Cj — 2ei/[ki(ki — 
1)]. Generally, in a network, for nodes with degree fc, their 
clustering coefficients, C(k), are not always the same. But 
in our network, all nodes with the same degree have identi- 
cal clustering coefficient. Moreover, for a single node with 
degree k, the analytical expression for its clustering coef- 
ficient C(k) can be derived exactly. 

According the connection rule (see figure [3]), when a 
node i enters the system, both ki and are 4. In the 
following steps, if one of its active triangles is selected, 
both ki and increase by 2 and 3, respectively. Thus, 
equals to 4 + | (ki — 4). The relation holds for all nodes 
at all steps. So one can see that there exists a one-to- 
one correspondence between the degree of a node and its 
clustering. For a node of degree k, we have 



C(k) 



2e 2[4+§(fc-4)] 



jfc(fc-l) 



k(k - 1) 



(4) 



In the limit of large k, C(k) exhibits a power-law behavior, 
C(k) ~ which has also been empirically observed in 
several real networks jTO] . 

We continue to compute the average clustering coeffi- 
cient C of RSN by means of the clustering spectrum C(k): 



C = J2P(k)C(k), 



(5) 
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Fig. 5. (Color online) Average clustering coefficient C of RSN 
vs the network size TV. The dotted line shows the analytic pre- 
diction and the squares denote the simulation results. 



6- 




N 

Fig. 6. (Color online) Semilogarithmic plot of the average path 
length L verse network size N. 



which can be easily obtained with respect to the degree 
distribution P{k) expressed by Eq. ([3]). The result is 



32 



k 

32 In 2 



k{k + 2)(k-2 
4tt2 

it 



1 



k- 1 k-1 



19 

— w 0.5674. 

* > 
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Thus the average clustering coefficient C of RAN is large 
and independent of network size. We have performed ex- 
tensive numerical simulations of the RSN. In figure OH we 
present the simulation results about the average cluster- 
ing coefficient of RSN, which are in complete agreement 
with the analytical value. 



4.3 Average path length 

From above discussions, we find that the existing model 
shows both the scale- free nature and the high clustering at 
the same time. In fact, our model also exhibits small-word 
property. Next, we will show that our network has at most 
a logarithmic average path length (APL) with the number 
of nodes. Here APL means the minimum number of edges 
connecting a pair of nodes, averaged over all couples of 
nodes. 

Using an mean-field approach similar to that presented 
in Ref. [20] , one can predict the APL of our network ana- 
lytically. By construction, at each time step, three nodes 
are added into the network. In order to distinguish differ- 
ent nodes, we construct a node sequence in the following 
way: when three new nodes are created at a given time 
step, we label them as M + 1, M + 2, . . . , M + 3, where 
M is the total number of the pre-existing nodes. Eventu- 
ally, every node is labeled by a unique integer, and the 
total number of nodes is N t — it + 3 at time t. We denote 
L(N) as the APL of our network with size N. It follows 

that L(N) = ^ N) 1} , where D{N) = T,i<i<j<N d i,j is 
the total distance, and where d%j is the smallest distance 



between node i and j. Note that the distances between 
existing node pairs are not affected by the addition of 
new nodes. As in the analysis of [5D], we can easily derive 
that D(N) ~ N 2 \nN in the infinite limit of N. Then, 
L(N) ~ InN. Thus, there is a slow growth of the APL 
with the network size N. This logarithmic scaling of L(N) 
with network size N, together with the large clustering co- 
efficient obtained in the preceding subsection, shows that 
the considered graph has a small- world effect. In figure [H 
we report average path length L(N) versus network size 
N. One can obviously see that L(N) increases logarithmi- 
cally with N. 



4.4 degree correlations 

First, we study the time evolution for the connectivity 
of an arbitrary node. Notice that the growing precess of 
RSN actually contains the preferential attachment mech- 
anism, which arises in it not because of some special rule 
including a function of degree as in Ref. [8] but natu- 
rally. Indeed, the probability that new nodes created at 
time t will be connected to an existing node i is clearly 
proportional to the number of active triangles contain- 



ing i, i.e. to its 



M*) 



Thus a node i is selected with 



the usual preferential attachment probability IIi[ki(t)] — 
[ki(t) - 2]/[2(2£ + 1)] - [h{t) - 2] /At (for large t). Conse- 
quently, ki satisfies the dynamical equation pQ: 



dkjjt) = h{t)-2 
dt At ' 

Considering the initial condition fcj(tj) = 4, we have 

ft\ 1/2 
fci(t) = 2(- +2. 



(7) 



(8) 



Having obtained the degrees for all nodes, we now 
study the degree correlations. Generally, degree correla- 
tions in a network can be conveniently measured by means 
of the quantity, called average nearest-neighbor degree (ANND), 
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which is a function of node degree, and is more convenient 
and practical in characterizing degree-degree correlations. 
The ANND is defined by [21] 



k nn {k) = ^2k'P(k'\k), 



(9) 



where P(k'\k) is the condition probability that a link be- 
longing to a node with connectivity k points to a node 
with connectivity kl . If there are no two degree correla- 
tions, k nn {k) is independent of k. When k nn (k) increases 
(or decreases) with k, the network is said to be assortative 
(or disassortative) [22] . 

Correlations can also be described by a Pearson corre- 
lation coefficient r, which is defined as l22l: 
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(10) 



where j m , k m are the degrees of the vertices at the ends 
of the mth edge, with m = 1, 2, • • ■ , M, where M denotes 
the number of edges in the network. The coefficient is in 
the range — 1 < r < 1. If the network is uncorrelated, the 
correlation coefficient equals zero. Disassortative networks 
have r < 0, while assortative graphs have a value of r > 0. 

We can analytically calculate the function value of 
k nn {k) for the RSN. Let Ri(t) denote the sum of the de- 
grees of the neighbors of node i, evaluated at time t. It is 
represented as 

jefi(i) 

where fl(i) corresponds to the set of neighbors of node 
i. The ANND of node i at time t, k nn (i,t), is then given 
by k nn (i,t) = Ri(t)/ki(t). During the growth of the RSN, 
Ri (t) can only increase by the addition of new nodes con- 
nected either directly to i, or to one of the neighbors of i. 
In the first case Ri(t) increases by 8 (the sum of degree 
for two newly-created nodes), while in the second case it 
increases by 2. Therefore, in the continuous k approxima- 
tion, we can write down the following rate equation [23j : 



dRj(t) 
dt 



8/7 i [fc i (t)] + 2 £ IljMt)] 

2[fc,(t)-2] Ri(t)-2ki(t) 
t 2t 

Ri(t) , h(t) 4 



2t t i 

The general solution of Eq. (|12[) is 

R l (t) = 4 + <P {i)ti +2 



1/2 



hit, 



(12) 



(13) 



where ^o(^) is determined by the boundary condition Ri(i). 
To obtain the boundary condition Ri(i), we observe that 
at time i, the new node i is connected to an existing node 



Fig. 7. (Color online) Plot of average nearest-neighbor degree 
of the nodes with degree k. The squares denote the simulation 
results, while the solid line is the theoretical result provided by 
Eq. (fT8]l. 



j of degree kj(i) with probability IIj[kj(i)], and that the 
degree of this node increase by 2 in the process. Thus, 



Ri(i) 



2] + 8, 



(14) 



where the last term 8 denotes the sum of the other two 
new nodes created at the same time as node i. Inserting 

/ . \ 1/2 

= 2 (jj +2 into Ri{i) 
leads to 



n 3 [kM = ^^F 1 and kj® 



Ri(i) = 8 + ln?: + E2(r/)~5 < 8 + 31ni. (15) 
i=i 

So, in the large i limit, Ri(i) is dominated by the second 
term, yielding 

Ri(i)<3lni. (16) 

From here, we have 



1/2 



Ri(t)~2 - lnt, (17) 



and finally 

k nn (k,t) ~lni. (18) 

So, two node correlations do not depend on the degree. 
The ANND grows with the network size N ~ 3i as In TV, 
in the same way as in the Barabasi- Albert (BA) model [24 
and the two-dimensional random Apollonian network |25j . 
In order to confirm the validity of the obtained analytical 
prediction of ANND, we performed extensive numerical 
simulations of the RSN (see figure [7]) with order N = 
100002. To reduce the effect of fluctuation on simulation 
results, the simulation results are average over fifty net- 
work realizations. From figure [7] we observe that for large 
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Fig. 8. (Color online) Semilogarithmic graph of Pearson cor- 
relation coefficient r as a function of network size N. 



k the ANND of numerical and analytical results are in 
agreement with each other, while the simulated results 
of ANND of small k have a very weak dependence on k, 
which is similar to the phenomena observed in the BA 
model [21] • This k dependence, for small degree, cannot 
be detected by rate equation approach, since it has been 
formulated in the continuous degree k approximation. 

To further confirm that RSN is uncorrelated, we com- 
pute the Pearson correlation coefficient r according to 
Eq. (|10p. The numerical results are reported in figure [5J 
From this figure, we see that for networks with small size, 
r is negative and only a little smaller than zero; when the 
size of the network increases, r goes to zero and is inde- 
pendent of size. The phenomenon of the convergence of r 
to zero again indicates that RSN shows absence of degree 
correlations. 

The fact that there is no degree correlation in RSN 
is compared with that of its deterministic variant, which 
shows a negative degree correlation [11]. We guess that 
the reason behind this difference between RSN and its de- 
terministic counterpart may stem from a biased choice of 
"active triangles" . In the evolution process of RSN, only 
one active triangle is updated at each iteration, while for 
the deterministic version, all active triangles are updated 
in one iteration. The difference between asynchronous and 
synchronous updating leads to different power-law expo- 
nents of degree distribution for the two networks, see [25] 
for detailed explanation. We think that the different as- 
sortative nature between the two graphs might be also 
related to a biased choice of active triangles at each it- 
eration. Of course, the genuine reason for this difference 
requires further study. 



4.5 modularity 

Many social and biological networks are fundamentally 
modular [T0l[27l[28] . These networks are formed by com- 
munities (modules) of nodes that are highly interconnected 
with each other, but have only a few or no links to nodes 
outside of the community to which they belong to. The 
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Fig. 9. (Color online) Modularity Q as a function of number 
of communities N q . r 
split with largest Q. 



strength of community structure is quantified through the 
modularity [29] 



8=1 



k 

E 



ds_ 
2E 



(19) 



where the sum runs over all communities, N q is the num- 
ber of communities (modules), E is the link number in the 
network, l s is the total number of links in the sth com- 
munity, and d s is the sum of the connectivities (degrees) 
of the nodes in module s. The modularity is high if the 
number of within- community links is much larger than 
expected from chance alone. 

We now look at the community structure of the net- 
work using the algorithm originally proposed by Girvan 
and Newman (GN) to find the partition with the largest 
modularity [lOj . The GN algorithm works by beginning 
with the complete network and at each step removing the 
edge with largest betweenness, where this quantity is re- 
calculated after the removal of every edge. If there is more 
than one edge with the same largest betweenness, we re- 
move them all at the same step. After all edges are re- 
moved, the network breaks up into N t communities (non- 
connected nodes). It is of interest to examine how the net- 
work is progressively broken into separate communities as 
one removes more edges. In figure [H we present how Q 
varies as the complete network with size N = 600 is bro- 
ken up into communities. Obviously, Q has a broad large 
value as a function of the number of communities, show- 
ing pronounced modular structure. Q is greater than 0.5 
when there are between 3 and 230 communities, and the 
best division has N q — 9 with Q = 0.7774. This largest 
value of Q is in contrast to the highest values found for 
some real-life networks such as collaboration network of 
scientists. Interestingly, the combination of modular and 
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uncorrelated properties has never been reported in previ- 
ous models. 



5 Conclusion 

In this paper, we have introduced a stochastic Sierpinski 
gasket and related it to a random maximal network, called 
random Sierpinski network (RSN). We have also proposed 
a iterative algorithm generating RSN, based on which we 
have determined some relevant topological characteristics 
of the network. We have presented that the network is 
simultaneously scale-free, small-world, uncorrelated, and 
modular. Thus, the RSN successfully reproduces some re- 
markable properties of many natural and man-made sys- 
tems. Our study provides a paradigm of representation for 
the complexity of many real-life systems, making it pos- 
sible to study the complexity of these systems within the 
framework of network theory. 

Acknowledgment 

We thank Yichao Zhang for preparing this manuscript. 
This research was supported by the National Basic Re- 
search Program of China under grant No. 2007CB310806, 
the National Natural Science Foundation of China under 
Grant Nos. 60496327, 60573183, 60773123, and 60704044, 
the Shanghai Natural Science Foundation under Grant 
No. 06ZR14013, the Postdoctoral Science Foundation of 
China under Grant No. 20060400162, Shanghai Leading 
Academic Discipline Project No. B114, the Program for 
New Century Excellent Talents in University of China 
(NCET-06-0376), and the Huawei Foundation of Science 
and Technology (YJCB2007031IN). 



References 

1. R. Albert and A.-L. Barabasi, Rev. Mod. Phys. 74, 47 

(2002) ; S. N. Dorogovtsev and J.F.F. Mendes, Adv. Phys. 
51, 1079 (2002); M.E.J. Newman, SIAM Rev. 45, 167 

(2003) ; S. Boccaletti, V. Latora, Y. Moreno, M. Chavez 
and D.-U. Hwanga, Phys. Rep. 424, 175 (2006). 

2. M. Faloutsos, P. Faloutsos, C. Faloutsos, Comput. Com- 
mun. Rev. 29, 251 (1999). 

3. R. Albert, H. Jeong, A.-L. Barabasi, Nature 401, 130 
(1999). 

4. H. Jeong, B. Tombor, R. Albert, Z.N. Oltvai, A.-L. 
Barabasi, Nature 407, 651 (2000). 

5. H. Jeong, S. Mason, A.-L. Barabasi, Z. N. Oltvai, Nature 
411, 41 (2001). 

6. M.E.J. Newman, Proc. Natl. Acad. Sci. USA 98, 404 
(2001). 

7. F. Liljeros, C.R. Edling, L.A.N. Amaral, H.E. Stanley, Y. 
Aberg, Nature 411, 907 (2001). 

8. A.-L. Barabasi and R. Albert, Science 286, 509 (1999). 

9. D.J. Watts and H. Strogatz, Nature (London) 393, 440 
(1998). 

10. M. Girvan and M. E. J. Newman, Proc. Natl. Acad. Sci. 
U.S.A. 99, 7821 (2002). 



11. Z. Z. Zhang, S. G. Zhou, T. Zou, L. C. Chen, and J. H. 
Guan, Eur. Phys. J. B 60, 259 (2007). 

12. W. Sierpinski, Comptes Rendus (Paris) 160, 302 (1915). 

13. C. A. Reiter, Comput. Graphics 18, 885 (1994) 

14. S. Hutchinson, Indiana Univ. Math. J. 30, 713 (1981). 

15. D.B. West, Introduction to Graph Theory (Prentice-Hall, 
Upper Saddle River, NJ, 2001). 

16. J.S. Andrade Jr., H.J. Herrmann, R.F.S. Andrade and 
L.R.da Silva, Phys. Rev. Lett. 94, 018702 (2005); Z.Z. 
Zhang, S.G. Zhou, L.J. Fang, J.H. Guan, Y.C. Zhang, EPL 
(Europhys. Lett.) 79 38007 (2007); Z.Z. Zhang, F. Cornel- 
las, G. Fertin, A. Raspaud, L.L. Rong, and S. G. Zhou, 
J. Phys. A: Math. Thero. 41, 035004 (2008); Z.Z. Zhang, 
L.C. Chen, S.G. Zhou, L.J. Fang, J.H. Guan, and T. Zou, 
Phys. Rev. E 77, 017102 (2008). 

17. PL. Krapivsky, S. Redner, and F. Leyvraz, Phys. Rev. 
Lett. 85, 4629 (2000). 

18. M. Hinczewski and A. N. Berker, Phys. Rev. E 73, 066126 
(2006); Z. Z. Zhang, S. G. Zhou, and T. Zou, Eur. Phys. 
J. B 56, 259 (2007); M. Hinczewski, Phys. Rev. E 75, 
061104 (2007); Z. Z. Zhang, S. G. Zhou, and L. C. Chen, 
Eur. Phys. J. B 58, 337 (2007); L. Wang, F. Du, H. P. Dai, 
and Y. X. Sun, Eur. Phys. J. B 53, 361 (2006). 

19. E. Ravasz and A.-L. Barabasi, Phys. Rev. E 67, 026112 

(2003) . 

20. T. Zhou, G. Yan, and B. H. Wang, Phys. Rev. E 71, 046141 
(2005); Z.Z. Zhang, L.L. Rong and F. Cornelias, Physica 
A 364, 610 (2006); L. Wang, H. P. Dai, and Y. X. Sun, J. 
Phys. A: Math. Thero. 40, 13279 (2007). 

21. R. Pastor-Satorras, A. Vazquez and A. Vespignani, Phys. 
Rev. Lett. 87, 258701 (2001). 

22. M. E. J. Newman, Phys. Rev. Lett. 89, 208701 (2002). 

23. A. Barrat and R. Pastor-Satorras, Phys. Rev. E 71, 036127 
(2005). 

24. A. Vazquez, R. Pastor-Satorras and A. Vespignani, Phys. 
Rev. E 65, 066130 (2002). 

25. Z. Z. Zhang and S. G. Zhou, Physica A 380, 621 (2007). 

26. F. Cornelias, H. D. Rozenfeld, D. ben-Avraham, Phys. Rev. 
E 72, 046142 (2005). 

27. E. Ravasz, A. L. Somera, D. A. Mongru, Z. N. Oltvai, and 
A.-L. Barabsi, Science 297, 1551 (2002). 

28. F. Radicchi, C. Castellano, F. Cecconi, V. Loreto, and D. 
Parisi, Proc. Natl. Acad. Sci. U.S.A. 101, 2658 (2004). 

29. M. E. J. Newman and M. Girvan, Phys. Rev. E 69, 026113 

(2004) . 



